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Based on the model of the Institute of Silicate Chemistry, a mathematical model of glass beam deformation is 
proposed. This model is tested using the experimental data, and its application makes it possible to predict the 
behavior of glass beam samples in molding. 


Multiple attempts to give an adequate mathematical de¬ 
scription of the sheet glass deformation process which takes 
place in production of bent articles often produced rather 
contradictory results [1 -3]. In every case there was a ques¬ 
tion with respect to the accuracy of the analytical calcula¬ 
tions which could be only verified experimentally. Such a 
possibility emerged when researchers at the Belgorod Tech¬ 
nological Academy of Construction Materials developed a 
viscometer which makes it possible to simulate cross-section 
flexural strain of a beam made from a particular glass com¬ 
position over a wide temperature range. 

The design of the instrument provided for strict obser¬ 
vance of the prescribed temperature conditions (with an ac¬ 
curacy of about 2%) and measurement of the current flexural 
strain of the beam using a mechanotron. Some results of the 
investigation of window pane plass were kindly presented to 
us by Yu. L. Belousov (Cand. Sci., Belgorod Technological 
Academy) and are reported in the present paper. 

A preliminary analysis of a number of mathematical 
models describing processes of temporary deformation of 
samples at high temperatures indicated that the model of the 
ISC (Instutute of Silicate Chemistry, Russian Academy of 
Sciences) developed precisely for glass [4] is the most ade¬ 
quate. However, its authors did not pose the problem of de¬ 
scribing the deformation of complex design schemes and 
limited themselves to the case of axial tension of the beam. 
In our study, we wanted to test the adequacy of the ISC 
model for the case of beam flexure under the effect of an in¬ 
tegrated load. 

In particular, the beam in Fig. 1 is loaded with distributed 
gravity load q and concentrated load F applied to the center 
of the span. This design scheme fully corresponds to the de¬ 


formation of the actual sample in the conditions of the exper¬ 
iment described above. 

The maximum value of the sag for purely elastic defor¬ 
mation will be [5] 
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where / is the beam length; E is the elongation modulus; J is 
the axial moment of inertia of the beam section. 

Since the ISC model describes the elongation variation 
with time x, i.e. 

s(x)=^—t + e (co)[l - M (t)] , 

n 


where a(0) is the initial strain at the time moment x = 0; r| is 
the material viscosity; £ (/ (co) is the creep component of de¬ 
formation at x -» oo; M £ (t) is the relaxation function 


A/ (x)= 


1-exp( 


-( T / T P ) 


where x is the current time; x r is the duration of stress relax¬ 
ation; b = 0.5, the relationship between the maximum sag 
(o max and the static elongation of the lower (stretched) fiber 
e was found: 
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Thereafter it is sufficiently easy to find the expression for 
determining the temporal variation in beam sag: 
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Thus, in order to calculate co(x), it is necessary first to de¬ 
termine r| = r|(r), £j(®) and x r = x r (f), where t is the temper¬ 
ature of glass. The viscous and relaxation parameters of glass 
for a given temperature level can be calculated according to 
O. B. Mazurin [6], and s^foo) can be found in the following 
way. 

Let us decompose plot e(x) into its components in the 
same coordinate system (Fig. 2). It can be seen that curve 
e(x) can be found from the geometric summation of e d (x) 
and £ n (x), where £ n (x) is the viscous deformation component 
calculated from the formula: 


With an increase in x at fixed temperature and strain val¬ 
ues, function £ rf (x) approaches a certain constant value 
£ (/ (x) m ax which is parameter s d (co), i.e., it can be stated that at 
a certain time, the further course of the function e(x) will de¬ 
pend only on £^(x), and the contribution of the creep £ d (x) to 


co(x) x 10\ mm 



Fig. 3. Variation in beam flexure. Solid curves) estimated values; 
dots) experimental data; 1) at temperature of 640°C, a(0) = 16.29 kPa; 
2) at temperature of 660°C, cr(0) = 17.09 kPa. 

the total deformation is restricted by the constant value 
£^( 00 ). A characteristic point in this case will be e(x ljm ), 
where for time x, im the values of the components £ n (x) and 
£ rf (x) will be equal. Let us find this point. 

Assuming that the value of z d (oo) for constant tempera¬ 
ture and loading conditions is constant, x hm can be calcu¬ 
lated, provided that the function 

f(t) = 1 - MJj) 

at this point has a known value. The value of function /(x) 
can be theoretically determined with sufficiently high preci¬ 
sion. Since the curve has substantial curvature and rapidly 
approaches an asymptote, it can be stated that the point x (jm 
will be positioned within the interval of the function with 
values above 0.95 when the further effect of £^(x) will be 
minimum. By denoting the accepted value of/(x) as C, it can 
be written: 


x„ m = x r (-ln(l-0) ,/ *. 

After finding x lim , it is possible to determine £ d ( co ) since 
for time x hm , the values of £ d (x) and £ n (x) will be equal (cal¬ 
culated in the same coordinate system): 

g(0)x p 
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Using the described method, the current flexure values of 
beams made of window pane glass were calculated for tern- 
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peratures of 540 and 660°C (Fig. 3). The sizes of the beams 
in the calculations corresponded to the values used in the ex¬ 
periments. At 640°C, / = 20 mm, b = 4.71 mm, h-3.12 mm; 
at 660°C, / = 24 mm, b = 4.678 mm, // = 3.415 mm. In addi¬ 
tion to that, the value of the initial stress o(0) of the beam 
was taken into account. 

The main criterion for verification of the model ade¬ 
quacy is the root-mean-square deviation calculated from the 
formula: 

i 

where 5 ; . is the deviation of the estimated values from the ex¬ 
perimental values; N is the number of measurements; k is the 
number of the degrees of freedom in the model. For the case 
considered, k = 3. 

As a result of the performed calculations, the following 
root-mean-square deviations were obtained: at the tempera¬ 
ture of 660°C, o = 17.09 kPa, ct 0 = 2,55; at the temperature 
of 640°C, a = 16.29 kPa, ct 0 = 2.176. 


As can be seen from the results obtained, the root-mean- 

square deviations are rather low, which is evidence of the ad¬ 
equacy of the proposed model to the actual deformation pro¬ 
cess in the glass beam. 
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